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Abstract 



J> ' In this paper, we calculate the branching ratios and CP-violating asymmetries for B° — » p°r]^ 

and B + — * p + n^ decays in the perturbative QCD factorization approach. In this approach, we 
not only calculate the usual factorizable contributions, but also evaluate the non-factorizable and 
annihilation type contributions. Besides the current-current operators, the contributions from 
', the QCD and electroweak penguin operators are also taken into account. The theoretical predic- 

tions for the branching ratios are Br(B + — > p + rj^) «9x 10~ 6 and Br(B° — > p°n^) ~ 5 x 10~ 8 , 
which agree well with the measured values and currently available experimental upper lim- 
its. We also predict large CP-violating asymmetries in these decays: A^jp(p ± n) ~ —13%, 
^Cp(pV) ~ -18%, ^cp(A) « -41%, A d j P (p°n') « -27%, A^(p°n) « +25%, and 
A?Q % p (p°n') ~ +11%, which can be tested by the current or future B factory experiments. 
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I. INTRODUCTION 



Along with the great progress in theoretical studies and experimental measurements, 
the charmless two-body B meson decays are getting more and more interesting and at- 
tracting more and more attentions, since they provide a good place for testing the standard 
model (SM), studying CP violation and searching for possible new physics beyond the 
SM. 

For the two-body hadronic B meson decays, the dominant theoretical error comes from 
the uncertainty in evaluating the hadronic matrix element < MiM 2 \Oi\B > where Mj 
and M2 are light final mesons. At present, the QCD factorization (QCDF) approach 0,01 
and the perturbative QCD (PQCD) factorization approach 0, 0, II are the two popular 
methods being used to calculate the hadronic matrix elements. The perturbative QCD 
approach has been developed earlier from the QCD hard-scattering approach Some 
elements of this approach are also present in the QCD factorization formula of Refs. 0, 0] • 
The two major differences between these two approaches are (a) the form factors are 
calculable perturbatively in PQCD approach, but taken as the input parameters extracted 
from other experimental measurements in the QCDF approach; and (b) the annihilation 
contributions are calculable and play an important role in producing CP violation for 
the considered decay modes in PQCD approach, but it could not be evaluated reliably 
in QCDF approach. Of course, one should remember that the assumptions behind the 
PQCD approach, specifically the possibility to calculate the form factors perturbatively, 
are still under discussion 0. More efforts are needed to clarify these problems. 

Up to now, many B meson hadronic decay channels have been calculated and studied 
phenomenologically in both the QCDF approach 0, 0, E| and in the PQCD approach 
|9l. lid. I 111 Il2l. Il3j|. In this paper, we would like to calculate the branching ratios and CP 
asymmetries for the B — ► pr/^ decays by employing the low energy effective Hamiltonian 
14 | and the PQCD approach. Besides the usual factorizable contributions, we here are 



able to evaluate the non-factorizable and the annihilation contributions to these decays. 

Theoretically, the four B — > prf'^ decays have been studied before in the naive or 
generalized factorization approach ^5|, as well as in the QCD factorization approach 
[8j]. On the experimental side, the branching ratios of B — > p + r],p + r]' decays have been 
measured 3, 17. lii lit]. 



Br(B + -> p + rj) = (8.lt\i) x 1(T 6 , 

Br(B + -> p+n') = (I2.9l|| ± 2.0) x 1(T 6 . (1) 



For B — > p°r], p°r]' decays, only the experimental upper limits are available now \li 

Br(B° -> p°r]) < 1.5 x 1(T 6 , Br(B° -> p°rj') < 4.3 x 1(T 6 . (2) 

In B — > pr]('> decays, the B meson is heavy, setting at rest and decaying into two light 
mesons (i.e. p and r/^ ) with large momenta. Therefore the light final state mesons are 
moving very fast in the rest frame of B meson. In this case, the short distance hard 
process dominates the decay amplitude. We assume that the soft final state interaction is 
not important for such decays, since there is not enough time for light mesons to exchange 
soft gluons. Therefore, it makes the PQCD reliable in calculating the B — > prf 1 ' decays. 
With the Sudakov resummation, we can include the leading double logarithms for all 
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loop diagrams, in association with the soft contribution. Unlike the usual factorization 
approach, the hard part of the PQCD approach consists of six quarks rather than four. 
We thus call it six-quark operators or six-quark effective theory. Applying the six-quark 
effective theory to B meson decays, we need meson wave functions for the hadronization of 
quarks into mesons. All the collinear dynamics are included in the meson wave functions. 

This paper is organized as follows. In Sec. |Hj we give a brief review for the PQCD fac- 
torization approach. In Sec. we calculate analytically the related Feynman diagrams 
and present the various decay amplitudes for the studied decay modes. In Sec. I1V[ we 
show the numerical results for the branching ratios and CP asymmetries of B — > pv^ > de- 
cays and comparing them with the measured values. The summary and some discussions 
are included in the final section. 



II. THEORETICAL FRAMEWORK 

The three scale PQCD factorization approach has been developed and applied in the 
non-leptonic B meson decays 0, 0, IE B EH U, El] for some time. In this approach, 
the decay amplitude is separated into soft, hard, and harder dynamics characterized by 
different energy scales (t, nib, iWjv). It is conceptually written as the convolution, 

A{B - M X M 2 ) ~ Jd 4 hd 4 k 2 d% Tr [C(t)^ B (k 1 )^ Ml (k 2 )^M 2 (k 3 )H(k 1 ,k 2 ,k 3 ,t)} , (3) 

where ki's are momenta of light quarks included in each mesons, and Tr denotes the trace 
over Dirac and color indices. C (t) is the Wilson coefficient which results from the radiative 
corrections at short distance. In the above convolution, C(t) includes the harder dynamics 
at larger scale than Mb scale and describes the evolution of local 4-Fermi operators from 
my/ (the W boson mass) down to t ~ 0{-\fAM B ) scale, where A = Mb— m^. The function 
H(ki, k 2 , k 3 ,t) describes the four quark operator and the spectator quark connected by a 
hard gluon whose q 2 is in the order of AM B , and includes the 0(\/AMb) hard dynamics. 
Therefore, this hard part H can be perturbatively calculated. The function $^ is the 
wave function which describes hadronization of the quark and anti-quark to the meson 
M. While the function H depends on the processes considered, the wave function $a/ is 
independent of the specific processes. Using the wave functions determined from other 
well measured processes, one can make quantitative predictions here. 

Since the b quark is rather heavy we consider the B meson at rest for simplicity. It is 
convenient to use light-cone coordinate {p + ,p~,Pr) to describe the meson's momenta, 

P ± = ^{P°±P 3 ), and p T = {p\p 2 ). (4) 

Through out this paper, we use the light-cone coordinates to write the four momentum as 
(kf, k± , k^). Using these coordinates the B meson and the two final state meson momenta 
can be written as 

Mr Mr n Mr 

A = ^f(M,0 T ), P 2 = -|(l,rJ,0 T ), P 3 = ^f(0,l,0 T ), (5) 

respectively, where r p = m p /mB\ the light pseudoscalar meson masses have been ne- 
glected. 
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For the B — ► prf^ decays considered here, only the p meson's longitudinal part con- 
tributes to the decays, its polar vector is €l = J^^ (1, — r 2 p , T )- Putting the light (anti-) 

quark momenta in B, p and 77 W mesons as fci, ki, and £3, respectively, we can choose 

h = (x 1 P+,0,k 1 T), k 2 = (x 2 P+,0,k 2T ), k 3 = (0,x 3 PzMt). (6) 
Then, the integration over k± , , and k 3 in eq.© will lead to 

A(B — > pr/'-*) ~ Jdxidx2dx 3 bidbib2db2b3db 3 

Tr [C(*)$ B (zi, p(^2, 63)^(0 ( x 3' & 3)^(^, 6*, t)^(x 4 ) e- 5 ^] ,(7) 

where 6j is the conjugate space coordinate of kiT, and t is the largest energy scale in 
function H(xi, 6j, £). The large logarithms (In ^? SL ) coming from QCD radiative corrections 
to four quark operators are included in the Wilson coefficients C{t). The large double 
loga rithms (In 2 Xi) on the longitudinal direction are summed by the threshold resummation 
[2(1 ] . and they lead to S t (xi) which smears the end-point singularities on x«. The last term, 
e -S(t)^ j g Sudakov form factor resulting from overlap of soft and collinear divergences, 
which suppresses the soft dynamics effectively [21]. Thus it makes the perturbative 
calculation of the hard part H applicable at intermediate scale, i.e., Mb scale. We will 
calculate analytically the function H(xi, bi, t) for B — > prp') decays in the first order in a s 
expansion and give the convoluted amplitudes in next section. 



Wilson coefficients 



For B — > pr/W decays, the related weak effective Hamiltonian H e ff can be written as 



n 



eff 



Gp 
71 
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v U bV* d (0^)0^) + c 2 (p)o u M) - v tb v; d £ c^) o^) 



i=3 



We specify below the operators in H. e ff for b — > d transition: 



O'l = dyflMp ■ upj^Lba , 0% 

3 = da^Lba ■ Eg' Qpl^Qp , 4 



doff Lbp ■ Eg' q'f3l» L <la , 

\d a ^Lbfi ■ Y ql e q/ q'^^Rq' a 
\d a YLbp ■ E 9 / eq'q'p^Lq'a 



(9) 



where a and /3 are the SU(3) color indices; L and i? are the left- and right-handed 
projection operators with L = (1 — 75), R = (1 + 75). The sum over q' runs over the 
quark fields that are active at the scale p = 0(m b ), i.e., (q'e{u, d, s, c, b}). The PQCD 
approach works well for the leading twist approximation and leading double logarithm 
summation. For the Wilson coefficients Cj(p) (i = 1, . . . , 10), we will also use the leading 
order (LO) expressions, although the next-to-leading order calculations already exist in 
the literature 0] • This is the consistent way to cancel the explicit p dependence in the 
theoretical formulae. 
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For the renormalization group evolution of the Wilson coefficients from higher scale to 
lower scale, we use the formulae as given in Ref. directly. At the high my/ scale, the 
leading order Wilson coefficients Ci(Mw) are simple and can be found easily in Ref. • 



In PQCD approach, the scale t is chosen at the maximum value of various subprocess 
scales to suppress the higher order corrections, which may be larger or smaller than the 
rrib scale. In the range of m& < t < my/, we will evaluate the Wilson coefficients Cj(i) at 
scale t by using the leading logarithm running equations, as given explicitly in Eq.(Cl) 

of Ref. 0. In numerical calculations, we also use a s = 47r/[/3i hx{t 2 / a!q CD )] which is 

the leading order expression with A { q ] cd = 193MeV, derived from A ( q ] cd = 250MeV. Here 
(3\ = (33 — 2nf)/12, with the appropriate number of active quarks nf. rif = 5 for t > m b . 

At a given energy scale t = nib = 4.8 GeV, the LO Wilson coefficients Ci{m b ) as given 
in Ref. [9| are 

Ci = -0.2703, C 2 = 1.1188, C 3 = 0.0126, C 4 = -0.0270, 

C 5 = 0.0085, C 6 = -0.0326, C 7 = 0.0011, C 8 = 0.0004, 

C 9 = -0.0090, Cio = 0.0022. (10) 

In the range of t < rrib, then we evaluate the Wilson coefficients Cj(t) by using the 
Ci(rrib) in Eq. (fTU|) as boundary input and the leading logarithmic running equations as 
given in Appendix D of Ref. [9J] for the case of rif = 4. For the Wilson coefficient C^i), 
for example, the running equation in the two different regions can be written as 

C 2 (t) = \ (tT 6 / 23 + rr 2 / 23 ) , for m b <t< m w , (11) 
C 2 (t) = I (r^- 6 / 23 + r/ 2 / 23 ) (r 6/25 + e V25 ) 

+\ (r 6/23 - V 2/23 ) (r 6/25 - e 2/25 ) , for t< m b , (12) 

where rj = a>s{t)/a>s{rnw) and ^ = a>s{t) / a>s{rnb). For the running equations of other 
Wilson coefficients one can see Appendix C and D of ref. 



B. Wave Functions 



In the resummation procedures, the B meson is treated as a heavy-light system. In 
general, the B meson light-cone matrix element can be decomposed as fK 22l 



1 d 4 z 

o W7 

i 



3 iki -z 



(01^(0)^(^15^5)) 



(ki) 



V2 



:ki) 



(13) 



8a 



where n = (1,0, Ot), and v = (0, 1, Ot) are the unit vectors pointing to the plus and 
minus directions, respectively. From the above equation, one can see that there are two 
Lorentz structures in the B meson distribution amplitudes. They obey to the following 
normalization conditions 



d 4 h 



>B 



(ki 



ft 



d 4 h - 



>B 



(ki 



0. 



(14) 
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In general, one should consider these two Lorentz structures in calculations of B meson 
decays. However, it can be argued that the contribution of 4>b is numerically small 
j23l ]. thus its contribution can be numerically neglected. Using this approximation, we 
can reduce one input parameter in our calculation. Therefore, we only consider the 
contribution of Lorentz structure 

$b = ->==(^b + m B ) 75 i j(k 1 ). (15) 

In the next section, we will see that the hard part is always independent of one of the 
and/or k±, if we make approximations shown in next section. The B meson wave function 
is then the function of variable k± (or kf) and k±, 

MK,ki) = f dk+<p(kf, ki,ki). (16) 

The wave function for dd components in 77 W meson are given as: 
1 



* V JP, C) = -7= jtfjx) + + Cmg-W - v ■ n)^Jx) (17) 



where P and x are the momentum and the momentum fraction of r] d j, respectively. We 
assumed here that the wave function of ry d( j is same as the 7r wave function. The parameter 
( is either +1 or —1 depending on the assignment of the momentum fraction x. 

In B — > pr/W decays, p meson is longitudinally polarized. We only consider its wave 



function in longitudinal polarization 23, 25], 

< p-(P, e L )\d a (z)u p (0)\0 >= / dxe™ p - z {{ [i4 p {x) + m p ^ p (x)} + m^(x)}(18) 

V 2^V C Jo 

The second term in above equation is the leading twist wave function (twist-2), while the 
first and third terms are sub-leading twist (twist-3) wave functions. 

The transverse momentum k x is usually conveniently converted to the b parameter by 
Fourier transformation. The initial conditions of leasing twist <f>i(x), i = B,p,r),rf, are of 
non-perturbative origin, satisfying the normalization 

r 1 1 

<p t (x,b = 0)dx = ^=f i , (19) 
with fi the meson decay constants. 



III. PERTURB ATIVE CALCULATIONS 

In the previous section we have discussed the wave functions and Wilson coefficients 
of the amplitude in eq.®. In this section, we will calculate the hard part H(t). This 
part involves the four quark operators and the necessary hard gluon connecting the four 
quark operator and the spectator quark. We will show the whole amplitude for each 
diagram including wave functions. Similar to the B — ► up decays 26\, the eight diagrams 
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contributing to the B — > pr]^ decays are shown in Figure 1. We first calculate the 
usual factorizable diagrams (a) and (b). Operators 0\, 2 , 3 , O4, O9, and 0\q are 
{V — A)(V — A) currents, the sum of their amplitudes is given as 

pi POO 

F ep = 4V2G F 7cC F m B / dx\dx 3 / bxdbib 3 db 3 (/>b(%i, 61) 

• { [(1 + x 3 )4> p (x 3 , 63) + (1 - 2s3)r / ,(^(x3, 63) + 4>fa 3 , b 3 ))] 
a a (tl)h e (xii X 3 , 61, 63) expf-^a^tg)] 

+2r p 0'(x 3 ,6 3 )a s (^)/ie(a:3 > a; 1 ,63,6i)exp[-5' a6 (^)]} , (20) 

where CV = 4/3 is a color factor. The explicit expressions of the functions h l e , the 
energy scales t\ and the Sudakov factors S a b(t) Can be found in the Appendix. In the 
above equation, we do not include the Wilson coefficients of the corresponding operators, 
which are process dependent. They will be shown later in this section for different decay 
channels. The diagrams Fig. 1(a) and 1(b) are also the diagrams for the B — > p form 
factor A Therefore we can extract A^ p from Eq. (JUJ). 

The operators O5, Oq, O7, and 0$ have a structure of (V — A)(V + A). In some decay 
channels, some of these operators contribute to the decay amplitude in a factorizable way. 
Since only the axial- vector part of (V + A) current contribute to the pseudo-sealer meson 
production, 

(p\V - A\B)(rj\V + A\0) = -(p\V - A\B)( V \V - A\0), (21) 

the result of these operators is opposite to Eq. (j2*Uj) . In some other cases, we need to do 
Fierz transformation for these operators to get right color structure for factorization to 
work. In this case, we get (S + P)(S — P) operators from (V — A) ( V + A) ones. For these 
(S + P)(S — P) operators, Fig. 1(a) and 1(b) give 

p 1 poo 

F e p p = 8V2G F irC F tfm B / dxidx 3 / bidbib 3 db 3 <pB{xi,bi) 

Jo Jo 

■ { [4> P (x 3 , b 3 ) + r p ((x 3 + 2)4> s p (x 3 , b 3 ) - s 3 0* (x 3 , 63))] 
■a s (t 1 e )h e (x 1 ,x 3 , h, 63) exp[-S ab (tl)] 

+ {xi<p p {x 3 , b 3 ) + 2r p (f) s p (x 3 ,b 3 )) a s (tl)h e (x 3 ,x 1 ,b 3 ,b 1 ) exp[-S ab (t 2 e )]} . (22) 

For the non-factorizable diagrams (c) and (d), all three meson wave functions are 
involved. The integration of b 3 can be performed using 5 function 5(b 3 — b 2 ), leaving only 
integration of b\ and b 2 . For the (V — A)(V — A) operators, the result is 

16 /*oo 

M ep = —G F nC F m B / dxidx 2 dx 3 / b 1 db 1 b 2 db 2 <I>b{xi, b 1 )^(x 2j b 2 ) 

v3 Jo Jo 

■ {x 3 [(j) p (x 3 ,b 2 ) -2r p (f) p (x 3 ,b 3 )] a s (t f )h f (x 1 , x 2 , x 3 ,b 1 ,b 2 ) exp[-S cd (t f )]} .(23) 

For the non-factorizable annihilation diagrams (e) and (f), again all three wave func- 
tions are involved. Here we have two kinds of contributions. M ao is the contribution 
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(a) 



B 




(b) 




(c) 



P B 




(d) 




,(') 



B 




B 




B 




,(') 



FIG. 1: Diagrams contributing to the B — ► pr/') decays (diagram (a) and (b) contribute to the 
-B — > p form factor A Q 



containing operator type (V — A)(V — A), while is the contribution containing oper 
ator type (V — A)(V + A). 

16 



M, 



ap 



—j=G ftiC / dx\dx 2 dx 3 / bidbib 2 db 2 (Pb(xi, 61) 
v3 io io 

■ { [a^pOca, b 2 )(j)^(x 2 , 6a) + r p r„ ((a; 3 - ar 2 ) (</^(>2, 62)^(^3, 62) + ^(^2, 62) 
■4>p(x 3 , 62)) + (xs + x 2 ) (<f>%(x2, b 2 )<f> p {x 3 , b 2 ) + <f%{x 2 , b 2 )<f) t p {x 3 , 62)))] 

■a s (t))h 1 f (xi,x 2 ,x 3 ,b 1 ,b 2 )exp[-S e f(t 1 f )) - [x 2 4> p (x 3 ,b 2 )(/)^(x 2 ,b 2 ) 

+r p r v ((x 2 - x 3 ) (<p*(x 2 , b 2 )(f) p (x 3l b 2 ) + <f%(x 2 , b 2 )(j) s p {x 3 , b 2 )) + r p r v 
■ ((2 + x 2 + x 3 )(j)^(x 2 , b 2 )(p s p (x 3 , b 2 ) -(2 — X2- x 3 )<f%(x 2 , b 2 )ct) p {x 3) 6 2 )))] 

■a s (t 2 f )h 2 f (x 1 ,x 2 ,x 3 ,b 1 ,b 2 )exp[-S ef (t 2 f )}} , (24) 
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where r v = r n = m^jmB- 
16 

Mq P = — j=G 'p7rC 'jrm^ / dx\dx 2 dx 3 / b±dbib 2 db 2 4>b{xi, b\) 
v3 Jo </o 

• { [a; 2 r„0p(x 3 , 6 2 ) (^(^2, ^2) + <f%(x 2 , b 2 )) - x 3 r p ((f) s p (x 3 , b 2 ) + 0*(x 3 , b 2 )) 

■<P^{x 2 , b 2 )}a s (t))h)(x 1 , x 2 , x 3 , b x , b 2 ) ex.p[-S e f(tf)) 

+ [(2 - x 2 )r v (j) p (x 3 , b 2 ) (rf(x 2 , b 2 ) + <^{x 2 , b 2 )) - (2 - x 3 )r p ((j) p (x 3 , b 2 ) 

+4> p (x 3 , b 2 )) (p^(x 2 , b 2 )] a s (t 2 f )h 2 f (xi, x 2 , x 3 , h, b 2 ) exp[-S ef (t 2 f )]} . (25) 

The f actor izable annihilation diagrams (g) and (h) involve only p and rf^ wave func- 
tions. There are also two kinds of decay amplitudes for these two diagrams. F ap is for 
{V — A)(V — A) type operators, and is for (S — P)(S + P) type operators, 

F ap = -4:V2~7i-G F C F fBm B / dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 

■ { [x 3 (j)p(x 3 , b 3 )(j)^(x 2 , b 2 ) + 2r p r v (j)^(x 2 , b 2 ) ((1 + x 3 )(j) p (x 3 , b 3 ) 
-(1 - x 3 )(f) p (x 3 , b 2 ))] a s (t? e )h a (x 2 , x 3 , b 2 , b 3 ) exp[-S gh (t 3 e )} 
+ [x 2 (p p (x 3 , b 3 )4>^(x 2 , b 2 ) + 2r p r 71 4> p (x 3 , b 3 ) ((1 + x 2 )4>^(x 2 , b 2 ) 
-(1 -x 2 )(f%(x 2 , b 2 ))] a s (tj)h a (x 3 ,x 2 ,b 3 ,b 2 )exp[-S gh (tj)]} , (26) 



/*1 poo 

F^p = -8V2G F nC F m B f B / dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 

■ { [2r v (f) p (x 3 , b 3 )(p^(x 2 , b 2 ) + x 3 r p (<f) s p {x 3 , b 3 ) - 0* (x 3 , b 2 )) <fr^(x 2 , b 2 )] 
■a s (t 3 e )h a (x 2 , x 3 , b 2 , 63) exp[-S gh (tl)) 

+ [2r p (p p (x 3 , b 3 )<f)^(x 2 , b 2 ) + x 2 r v (<fr^(x 2 , b 2 ) - 0j(x 2 , b 2 )) <fr p (x 3 , b 3 )] 
■a s (tl)h a (x 3 ,x 2 ,b 3 ,b 2 )exp[-S gh (tl))} . (27) 

In the above equations, we have assumed that x\ « x 2 ,x 3 . Since the light quark mo- 
mentum fraction x\ in B meson is peaked at the small x\ region, while quark momentum 
fraction x 2 of rf 1 ' is peaked around 0.5, this is not a bad approximation. The numerical 
results also show that this approximation makes very little difference in the final result. 
After using this approximation, all the diagrams are functions of = a^m^/v^ of B 
meson only, independent of the variable of kf. Therefore the integration of eq. (fT6*j) is 
performed safely. 

If we exchange the p and rf'' in Figure 1, the result will be different. Because this will 
switch the dominant contribution from B — > p form factor to B — > rf^ form factors. The 
new diagrams are shown in Figure 2. 

We firstly consider the factorizable diagrams Fig. 2(a) and 2(b). The decay amplitude 
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p p 




(g) (h) 



FIG. 2: Diagrams contributing to the B — > /j?/'** decays (diagram (a) and (b) contribute to the 
B — > 77W form factor Fq~*^ '). 

F e induced by inserting the (V — A)(V — A) operators is 

pi POO 

F e = 4V27rG F C F f p m 4 B / dxidx 3 / bidbib s db 3 <j)B{xi,bi) 

Jo Jo 

• { [(1 + x 3 )^(x 3 , h) + r„(l - 2x 3 ) «(x 3 , 63) + 0j> 3 , & 3 ))] 
•a s (tg)/i e (xi, x 3 , 61, 63) expf-S^tg)] 

+2r^(a;3,63)a s (^)^(^3,^,fe3,fci)exp[-S' o6 (^)]} . (28) 

These two diagrams are also responsible for the calculation of B — ► form factors F Q B ^ V 
and Fq^ v , These two form factors can be extracted from Eq. (|2*8j) . 

Since only the vector part of the (V + A) current contribute to the vector meson 
production, the decay amplitude F e p induced by inserting (V — A)(V + A) operators is 
identical with the amplitude F e as given in Eq. (J28|) . i.e., 

F e p = F e . (29) 
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Because neither scaler nor pseudo-sealer density gives contribution to a vector meson 
production, (p\S + P\0) = 0, we get F e 5+P = 0. 

For the non-factorizable diagrams Fig. 2(c) and 2(d), the corresponding decay ampli- 
tudes are 

16 f 1 f°° 

M e = -=G ' pi\C 'f^b / dx\dx 2 dx 3 / b\db^b 2 db 2 <Pb(xi, bi)(p p (x 2 , b 2 ) 

v3 Jo Jo 

■ {x 3 [(p^(x 3 ,b 2 ) - 2r v <f)T(x 3 , b 2 )] a 8 (t f )h f (xi, x 2 , x 3 , b u b 2 ) exp[-S cd (t f )}} ,(30) 



32 f 1 f°° 

M[ = —GFirCFr p m B / dx\dx 2 dx 3 \ b\db\b 2 db 2 $B{x\,b\) 

V 3 io Jo 

■ { [x 2 tf(x 3 , b 2 ) (<f) s p (x 2 , b 2 ) - 4> P (x 2 , b 2 )) + r v ((x 2 + x 3 ) (tf(x 3 , b 2 ) 
■<Pp(x 2 , b 2 ) + <f%{x 3 , b 2 )(j) p (x 2 , b 2 )) + (x 3 - x 2 ) ((j)^(x 3 , b 2 )(j) p (x 2 , b 2 ) 
+(frj l (x 3 ,b 2 )<f) s p (x 2 ,b 2 )))] a s (tf)hf(x 1 ,x 2 ,x 3 ,b 1 ,b 2 )exp[-S cd (t f )]} . (31) 

From the non-factorizable annihilation diagrams Fig. 2(e) and 2(f), we find the decay 
amplitude M a for (V - A)(V - A) operators, Mf for (V - A)(V + A) operators, 

16 /*oo 
M a = —=TrGFCFm B I dx\dx 2 dx 3 I bidbib 2 db 2 (j)B(xi,bi) 
v3 Jo Jo 

■ { [x 3 <f>p(x 2 , b 2 )(f)^(x 3 , b 2 ) + r p r v ((x 3 - x 2 ) ((p^(x 3 , b 2 )4> P {x 2 , b 2 ) + ^(x 3 , b 2 ) 
■<l>p{x2, b 2 )) + {x 2 + x 3 ) (0^(x 3 , b 2 )(f) p (x 2 , b 2 ) + 0^(x 3 , b 2 )(p P (x 2 , b 2 )))] 
•a s (t})/i}(xi, x 2 , x 3 , 6i, 6 2 ) exp[—S e f(if)] 

+ [x 2 (f> p (x 2 , b 2 )(p^(x 3 , b 2 ) + r p r v {(x 2 - x 3 ) (<j>%{x 3 , b 2 )4> P (x 2 , b 2 ) + 0^(x 3 , b 2 ) 
■<Pp(x 2 , b 2 )) + (2 + x 2 + x 3 )(/)P(x 3 , b 2 )(f) s p (x 2 , b 2 ) - (2 - x 2 - x 3 )^(x 3 , b 2 ) 
■(f) i p{x 2 ,b 2 ))]a s (t 2 f )h 2 f (x 1 ,x 2 ,x 3 ,b 1 ,b 2 )exp[-S ef (t 2 f )]} , (32) 



M; = M p ap . (33) 
For the factorizable annihilation diagrams Fig. 2(g) and 2(h), we have 

F a = -F ap , and F[ = -F[ p . (34) 

Now we are able to calculate perturbatively the form factors F ( f~ >,?( ) (0), Aq^ p (0), 
and the decay amplitudes for the Feynman diagrams after the integration over Xi and fej. 
When doing the above integrations over Xi and bi, we should include the corresponding 
Wilson coefficients Ci(tj) calculated at the appropriate scale tj using Eqs. (CI) and (Dl) 
of Ref . . Since we here calculated the form factors and amplitudes at the leading order 
( one order of a s (t)), the radiative corrections at the next order would emerge in terms of 
a s (t) ln(m/i), where m's denote some scales, like m#, 1/bi, . . ., in the hard part H(t). We 
select the largest energy scale among m's appearing in each diagram as the hard scale t's 
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for the purpose of at least killing the large logarithmic corrections partially, 



t\ = a t ■ max( v /xim jB , l/b h l/b 3 ) , 

t\ = a t - m&x(y/¥[m B , 1/6 X , l/b 3 ) , 

t\ = at - rrmx( y /x 3 ~m B , l/b 2 , 1/63) , 

t\ = a t - max(^m B , l/b 2 , 1/63) , 

t f = a t ■ max( v /xix 3 m s , y/x 2 x 3 m B , 1/bi, l/b 2 ) , 

t l f = a t ■ max( v /x 2 x 3 m B , l/6 2 ) , 

= a t • max( v / a;i + a; 2 + x 3 - - x 2 x 3 m B , y/x 2 x 3 m B , l/&i, I/&2) , (35) 

where the constant a t = 1 ± 0.1 is introduced in order to estimate the scale dependence 
of the theoretical predictions for the observables. 

Before we put the things together to write down the decay amplitudes for the studied 
decay modes, we give a brief discussion about the rj — rj' mixing and the gluonic component 
of the rf meson. 

The rj and r( are neutral pseudoscalar (J p = CT) mesons, and usually considered as 
mixtures of the SU(3)f singlet 7/1 and the octet rjs' 



with 



i] \ = f cosOp - sin 9 P \ ( r)$ 
rj J \ sin 8 P cos 9 P J \Vi 



■)]$ = —= (uu + dd — 2ss) 
V6 



(36) 



771 = —j= (uu + dd + ss) , (37) 
v3 



where 9 P is the mixing angle to be determined by various related experiments |2j|. From 
previous studies, one obtains the mixing angle 6 P between —20° to —10°. One best fit 
result as given in Ref. [28] is -17° < 9 P < -10°. 

As shown in Eqs. fl36|37j) . 77 and 77' are generally considered as a linear combination of 
light quark pairs. But it should be noted that the 77' meson may has a gluonic component 
in order to interpret the anomalously large branching ratios of B — ■> ^7/ and J/^f — > r] >r y 
2M 2^1 . In Refs. 28[ 2^, 3^1, the physical states 77 and 77' were defined as 



1 77 > = X v 
\rf > = Xrf 



uu + dd 



V2_ 
uu + dd 



V2 



+ Y v \ss >, 

+ Y v '\ss > +Z v >\gluonium >, 



(3* 



where X (,) , Y (,) and Z v > parameters describe the ratios of uu + dd, ss and gluonium 
(SU(3)f singlet) component of 77W, respectively. In Ref. [28], the author shows that the 
gluonic admixture in 77' can be as large as 26%, i.e. 

Z v ,/ (X v w + Y v(l) + Z n ,) < 0.26. (39) 
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According to paper |22|, a large SU(3) singlet contribution can help us to explain the 
large branching ratio for B — ► Krf decay, but also result in a large branching ratio for 
B -► K°7] decay, Br(B -> fT ^) ~ 7.0(13) x 10" 6 for 9 P = -20°(-10°) as given in Table 
II of Ref. ji^], which is clearly too large than currently available upper limits (l^ |: 



Br(B -> i^°r]) < 1.9 x 10" 



(40) 



Although a lot of studies have been done along this direction, but we currently still do 
not understand the anomalous gg — rf coupling clearly, and do not know how to calculate 
reliably the contributions induced by the gluonic component of rf meson. In this paper, 
we firstly assume that rf does not have the gluonic component, and set the quark content 
of r\ and rf as described by Eqs. f|36l37|) . We will also discuss the effects of a non-zero 
gluonic admixture of 7]' in next section. 

Combining the contributions from different diagrams, the total decay amplitude for 
B + — > p + rj decay can be written as 



£u ( Ci + 3^2 



6 



1 3 15 

~ 3 ~ 4 ~ 2 ~~ 2 3 



+C 



io 



f;Fi(o P ) - < 

7-C5 + C*6 — 77C7 

3 o 



TrCV + -C*8 — -Cg — -C10 



6 



C 8 Fi( 



6 



+M, 



-6 a 



i u c2 — it 



C% — -Cg 



C 3 + 2C 4 + 2C e 
1 



:C S 



+ -C 10 



-c 



10 



+(M ap + M e + M a ) [£ u Ci - &(C 3 + C 9 )] ■ F 1 (9 P ) 
-(2M a p p + Mf) Ct(C B + C r )-F 1 (e p ) 



+F e 



in ~Cl + C2 - 6 77^3 + C 4 + -C 9 + C 



10 



F±( 



Fi(9 p )\, (41) 



where £ u = V* b V ud , & = V t * b V td , and F^[9 P ) = - sin9 p + cos B p /y/2 and F 2 (0 P ) = - sin# p - 
v/2 cos^p are the mixing factors. The Wilson coefficients Ci should be calculated at the 
appropriate scale t using equations as given in the Appendices of Ref. |~ 
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Similarly, the decay amplitude for B — > p i] can be written as 



V6M(p r]) = F e 



1 3 
— C"\ — Ca ~\~ — Cj 
3 2 



ep 



e« ( Ci + -c 2 



1 5 

2^8 + + Cio 



Fx(e p ) 



1 



1 



+^ep £t ( 77^5 + C 6 - -CV - -Cg ) • ^1 



6 



—M, 



ep 



1 



^uC*2 — it ( C3 + 2C4 + 2C*6 + -C*8 — 2^ 9 2^* 10 



( C 4 + C 6 - -C 8 - -Cio j l- ] 2\"i>) 



3 13 
iuC 2 — it ( — C3 + 2^8 + 2^ 10 



+ (M ap + M a ) 
- (Mf + 2Mf ) 6 f C 5 - ±C 7 ) F 1 (6 P 



3 13 
iuC 2 — C,t [ —C3 — -C 8 + -Cg + 2^10 



(42) 



The decay amplitudes for B — > p + r/ and S — > p ?/ can be obtained easily from Eqs. (j41j) 
and (J4*2*j) by the following replacements 



J T]1 J 7} 

Fi(9 p ) 

F2(e P ) 



ed fs 



F[( 



cos + 



sin^p 

~7T 



F^e p ) 



cos sin 



(43) 



Note that the possible gluonic component of 7/ meson has been neglected here. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 



A. Input parameters and wave functions 

We use the following input parameters in the numerical calculations 



A. 



MS 



m. 



Vdd 



250MeV, U = 130MeV, f B = 190MeV, 
1.4GeV, f p = 200MeV, f K = 160MeV, 



Mf 



5.2792GeV, M, 



w 



80.41GeV. 



(44) 
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The central values of the CKM matrix elements to be used in numerical calculations are 

El 



(45) 



\v ud \ 


= 0.9745, 


\v ub \ 


= 0.0040, 


\Vtb\ 


= 0.9990, 


\Vtd\ 


= 0.0075. 



For the B meson wave function, we adopt the model 
(j>B(x,b) = N B x 2 (l — x) 2 exp 



M 2 B x 2 



2ul 



(46) 



where uj\, is a free parameter and we take Ub = 0.4 ± 0.04 GeV in numerical calculations, 
and Nb = 91.745 is the normalization factor for uj^ = 0.4. This is the same wave functions 
as in Refs. 0, EU H3 , which is a best fit for most of the measured hadronic B decays. 

For the light meson wave function, we neglect the b dependant part, which is not 
important in numerical analysis. We choose the wave function of p meson similar to the 
pion case [25t | 

3 



(j) p (x) 



^fpx(l - x) 1 + 0.18C 2 3/2 (2a; - 1) 
V6 



{ p 

2V6 



{3(2x - l) 2 + 0.3(2x - l) 2 [5 (2a; - l) 2 - 3] 



+0.21[3 - 30(2x - l) 2 + 35(2s - l) 4 ]} 



<f>' p {x) = ^/p T (l - 2x) [1 + 0.76(10x 2 - lOx + 1)] . 
The Gegenbauer polynomial is defined by 

Cl'\t) = I (5t 2 - 1) . 



(47) 

(48) 
(49) 

(50) 



For 77 meson's wave function, (j>^ and (fi^ represent the axial vector, pseudoscalar 
and tensor components of the wave function respectively, for which we utilize the result 



from the light-cone sum rule |3J| including twist-3 contribution: 
3 



00 



+a^y [21(1 - 2x) 4 - 14(1 - 2x) 2 + l] J , 

) [3(1 - 2x) 2 - 1] 



2^/, {l + ^(30^3 -\ P % di 



1 

+ - 



27 2 81 2 vm. 
-Smuj-i p„ p a, 

20 dd 10 dd 



J [35(1 - 2x) 4 - 30(1 - 2x) 2 + 3] I 



x 



V2N< 
1 



m - 2x) 



g + ( 5r ?3 - 



^-Ip^)(Wx 2 -10x + 1) 



(51) 
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with 



^Vdd 

Pv dS 



0.44, a 

m n /m l 



Vdd 
4 

Vdd 

) 



= 0.25, a 
m = 0.015, 



K 



4 



= 0.20, 
u 3 = -3.0. 



0.25, 



(52) 



We assume that the wave function of uu is same as the wave function of dd. For the 
wave function of the ss components, we also use the same form as dd but with m s s and f v 
instead of and f x , respectively. For f x and f y , we use the values as given in Ref. [22j 
where isospin symmetry is assumed for f x and SU(3) breaking effect is included for f y : 



fx fn j fy 



VI - fl 



(53) 



These values are translated to the values in the two mixing angle method, which is 
often used in vacuum saturation approach as: 



U = 169MeV, fi = 151MeV, 

# 8 = -25.9°(-18.9°), 0i = -7.1°(-0.1°), 



(54) 



where the pseudoscalar mixing angle 8 P is taken as —17° (—10°) |28|. The parameters 
(* = Vdd(uu),Vss) are defined as: 



m. 



^ldd(uu) 

o 



mz 



= m' = 



[m u + m d ) ' 



ml 3 " = 



2M 2 K 



mz 



(2m s 



(55) 



We include full expression of twist— 3 wave functions for light mesons. The twist— 3 
wave functions are also adopted from QCD sum rule calculations j33|. We will see later 
that this set of parameters will give good results for B — > pr/^ decays. Using the above 
chosen wave functions and the central values of relevant input parameters, we find the nu- 
merical values of the corresponding form factors at zero momentum transfer from Eqs.()20j) 
and 



4 s " 




= 0) 


= 0.37, 


F B - 




= 0) 


= 0.15, 


p B- 
^0 


v (g 2 


= 0) 


= 0.14. 



(56) 



These values agree well with those as given in Refs. [31|, |32j, |3 



B. Branching ratios 

For B — > pr]('' decays, the decay amplitudes in Eqs. (}4~Tj) and can be rewritten as 
M = V* b V ud T - V t * b V td P = V* b V ud T [1 + zS ***] , (57) 

where 

v t tv td 



v* b v ud 



(5f 
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is the ratio of penguin to tree contributions, a = arg — y* d y» is the weak phase (one of 

the three CKM angles), and 5 is the relative strong phase between tree (T) and penguin (P) 
diagrams. The ratio z and the strong phase 5 can be calculated in our PQCD approach. 
One can leave the CKM angle a as a free parameter and explore the CP asymmetry 
parameter dependence on it. 

For B — > p + f] decay, for example, one can find "T" and "P" terms by comparing the 
decay amplitude as defined in Eq. (}4"Tj) with that in Eq. (|5T|). 

T( P + V) = ^ " [F ep (d + ^d) + M ep C 2 

+F e Qd + d) + (M a + M e + M ap ) d | , (59) 
P(P + V) = ^ ■ {Fe P (-\c* - C 4 - \c 7 - ±C a + 5 -C 9 + do) ft 



F ep [ + C 6 - -C 7 - -Cg 



+M ep f-C 3 + 2C 6 - ^C 8 + ^C 9 + ^C 10 

+F e Qc 3 + C A + l -C, + do 

+ (M ap + M e + M a ) (d + d) + (2M a p + Mf) (d + d)} 
F 2 (6 n ) f„ /1„ 1„ 1„ 1 



V3 



{^ep Qd + \c 8 - \c, - ido) • /, s 



+M ep (d + d - \c % - ido) } • (60) 

Similarly, one can obtain the expressions of the corresponding tree and penguin terms for 
the remaining three decays. 

Using the "T" and "P" terms, one can calculate the ratio z and the strong phase 5 for 
the decay in study. For B + — > p + rj and p + rj' decays, we find numerically that 

z(p + rj) = 0.10, 5(p + r]) = -137°, (61) 
z(p + rf) = 0.15, 5(pV) = -139°. (62) 

The errors of the ratio z and the strong phase 5 induced by the uncertainty of the input 
parameters, such as u b = 0.4 ± 0.04 GeV, = 1.4 ± 0.1 GeV, and a = 100° ± 20°, are 
very small in magnitude and not be shown explicitly in Eqs. ()61|) and (|62|). The reason is 
that the errors induced by the uncertainties of these input parameters are canceled almost 
completely in the ratio. 

Unlike the case of QCD factorization approach, the energy scale t (in PQCD factor- 
ization approach ) appeared in the Wilson coefficients C»(0 an d in the Sudakov form 
factors Sj(t) vary simultaneously during the integration over Xj and &j (i = 1,2,3). If 
we choose the hard energy scale t'.s as defined in Eqs. (J35|) with a t = 1, there will be no 
remaining scale dependence left explicitly after the integration. But we know that such 
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scale dependence should exist and likely dominate the errors on theoretical predictions for 
those observables. Since the calculation in this paper is performed at the leading order 
and thus may suffers from the uncertainties due to the scale dependence of the LO Wilson 
coefficients. In Ref. (37j the authors calculated the branching ratios of B — > Kix, tttt firstly 
at the next-to-leading order by using the PQCD factorization approach, and they found 
that the NLO contribution can give about 15 — 20% correction to LO predictions. The 
size of NLO contribution in PQCD approach is indeed very complicated to calculate. To 
explore it, as shown in Eq. (|35|). we here mult iply a factor a t = 1 ± 0.1 to the ordinary 
definition of scale t,'s in Refs. U EH 11, EH 13l |. and take it as an estimation for the 



uncertainty of the possible scale dependence. Numerically we find that 

z(p + rj) = 0.10j£» 5(p + r ] ) = (-137i»)°, (63) 
^V) = 0.15±™i, = (-139^ 7 )°, (64) 

for a t = 1 ± 0.1. The larger change of z and 5 corresponds to the case of a t = 0.9, while 
the magnitude of the variations is consistent with the general expectation. 

From Eq. (JSTj) . it is easy to write the decay amplitude for the corresponding charge 
conjugated decay mode 

M = V ub V* d T - V tb V td P = V ub V: d T [1 + ze^ a+ V] . (65) 

Therefore the CP-averaged branching ratio for B° — > pr/^ is 

Br = (\M\ 2 + \M\ 2 )/2 = \V ub V* d T\ 2 [1 + 2z cos a cos 5 + z 2 ] , (66) 

where the ratio z and the strong phase 5 have been defined in Eqs. (j57j) and (J58)) . It is 
easy to see that the CP-averaged branching ratio is a function of cos a for the given ratio 
z and the strong phase 5. This gives a potential method to determine the CKM angle a 
by measuring only the CP-averaged branching ratios with PQCD calculations. But one 
should know that the uncertainty of theory is so large as to make it unrealistic. 

Using the wave functions and the input parameters as specified in previous sections, it 
is straightforward to calculate the branching ratios for the four considered decays. The 
theoretical predictions in the PQCD approach for the branching ratios of the decays under 
consideration are the following 

Br { B + - p+ V ) = [8.5ljjS(u*)lgS(mS) ± 0.4(a)^j(a t )] x 10" 6 , (67) 

Br{ B + -> p+r/) = [8.7^;°(^)^; 9 7 (mJ)^; 5 7 (a)^i(a,)] x 10" 6 , (68) 

Br( B° - p\) = [0.024l°;^K)lX«) ± 0.002(a)l^(a t )] x 10~ 6 , (69) 

Br( B° - (Prf) = [0.061^^K)^;^K) ± 0.003(a)^(a t )] x lO" 6 , (70) 

for 6 P = -10°; and 

Br( B+ - p+rj) = [l0.6i^(a; 6 )^;°K) ± 0.5(a)^(a t )] x lO' 6 , (71) 

Br( B + -> p+n') = [6.5t?:i(w6) ± 0.6(m5) ± 0.5(a)+^(a t )] x lO" 6 , (72) 

Br(B°-A) = [0.042lHf 2 (a; 6 )±0.005(mJ)tHo4(«)^?2(«t)])><10- 6 , (73) 

Br( B° AO = [0.047i°;^(a; 6 )lHo6«) ± 0.001(c)l^™(a t )] x lO" 6 , (74) 
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FIG. 3: The a dependence of the branching ratios (in unit of 10 6 ) of B + — ► p + r/^ decay for 



1.4 GeV, 9 P 



-10°, u h = 0.36 GeV (dotted curve), 0.40 GeV(solid curve) and 0.44 



GeV (short-dashed curve). The gray band show the data. 



for r . 



-17°. The major errors are induced by the uncertainty of hard energy scale 



t, u b = 0.4 ± 0.04 GeV, ml = 1.4 ± 0.1 GeV and a = 100° ± 20°, respectively. It is 
easy to see that (a) the errors of the branching ratios induced by varying a t in the range 
of at = [0.9, 1.1] are less than 20% for the tree-dominated B — > p + r/^ decays; but can 
be significant for the penguin-dominated B — > p 0r q^ decays; and (b) the variations with 
respect to the central values are large (small) for the case of a t = 0.9 (a t = 1.1). This 
feature agrees with general expectations: when the scale t become smaller, the reliability 
of the perturbative calculation of the form factors in PQCD approach will become weak! 

It is easy to see that the PQCD predictions for the branching ratios of considered 
decays agree very well with the measured values or the upper limits as shown in Eqs.(fT|) 
and (0). For the four B — > pr/^ decays, the theoretical predictions for the CP-averaged 
branching ratios in the PQCD approach are well consistent with those given in the QCD 
factorization approach 



Br( 


B+- 


-* P +r l) = 


(9.41*1) x 


10~ 6 , 


Br( 


B + - 


> pV) = 




10~ 6 , 


Br 


(B°- 


- A) = 


(0.03r°;l 7 ) 


x 10~ 6 


Br( 


B°- 


- A') = 


(o.oii^) 


x 10~ 6 



(75) 
(76) 

where the individual errors have been added in quadrature. 

It is worth stressing that the theoretical predictions in the PQCD approach have large 
theoretical errors induced by our ignorance of NLO contributions, and the still large 
uncertainties of many input parameters. In our analysis, we consider the constraints on 
these parameters from analysis of other well measured decay channels. For example, the 
constraint l.lGeV < < 1.9GeV was obtained from the phenomenological studies for 
B — > 7T7T decays 9], while the constraint of a ~ 100° ± 20° were obtained by direct 
measurements or from the global fit [3, |3^, |36| . From numerical calculations, we get to 
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FIG. 4: The a dependence of the branching ratios (in unit of 10 6 ) of B + — > p + rj^ decays for 
io b = 0.4 GeV, 8 p = -10°, ml = 1.3 GeV (dotted curve), 1.4 GeV (solid curve) and 1.5 GeV 
(short-dashed curve). The gray band shows the data. 
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FIG. 5: The a dependence of the branching ratios (in unit of 10 6 ) of p°7] (solid curve) and p°rj' 
(dotted curve) decays for ml = 1.4 GeV,6> p = -10°, uj b = 0.40 GeV . 



know that the main errors come from the uncertainty of u b , mj, a, 8 P and the next-to- 
leading order contributions. 

In Figs.|21and|3J we present, respectively, the PQCD predictions of the branching ratios 
of B -> p + T] and p+rf decays for 6 P = 10°, u b = 0.4 ± 0.04 GeV, m£ = 1.4 ± 0.1 GeV and 
Fig. shows the a-dependence of the PQCD predictions of the branching 
10°, u b = 0.4 GeV, ml = 1.4 GeV and a = [0°, 180°]. 



a = [0°, 180° 
ratios of B - 



From the numerical results and the figures we observe that the PQCD predictions are 
very sensitive to the variations of u b and mj. The parameter ml originates from the 
chiral perturbation theory and have a value near 1 GeV. The ml parameter characterizes 
the relative size of twist 3 contribution to twist 2 contribution. Because of the chiral 
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enhancement of mj, the twist 3 contribution become comparable in size with the twist 2 
contribution. The branching ratios of Br(B — > pq^) are also sensitive to the parameter 
mj, but not as strong as the ui^ dependence. 



C. CP- violating asymmetries 

Now we turn to the evaluations of the CP-violating asymmetries of B — > pq^ decays 
in PQCD approach. For B + — > p + r] and B + — > p + rj' decays, the direct CP-violating 
asymmetries A C p can be defined as: 

A dir = Igl 2 - l^l 2 = 2z sin a sin 5 
CP \M\ 2 + \M\ 2 1 + 2zcosacos5 + z 2 ' 1 ' 

where the ratio z and the strong phase 5 have been defined in previous subsection and 
are calculable in PQCD approach. 

Using the central values of z and 6 as given in Eqs. lJoTj) and (jo^j) . it is easy to calculate 
the CP-violating asymmetries. In Fig. |Hl we show the a— dependence of the direct CP- 
violating asymmetries A d Q P for B ± — > p ± rj (the solid curve) and B ± — > p ± rj' (the dotted 
curve) decay, respectively. From Fig. El one can see that the CP- violating asymmetries 
A d Qp[B ± — >■ p ± T]) and Aq P {B ± — > p ±r q') are large in magnitude, about —15% for a ~ 100°. 
The large CP-violating asymmetries plus large branching ratios are clearly measurable in 
the B factory experiments. 

For a = 100° ± 20°, one can read out the allowed ranges of Afj p from Fig. Indirectly 

A&{B± ^ fa) = (-IStlH^tlM)) x 10- 2 , 

A&{B± - pS) = (-lZtH(a)±Ua t )) x (78) 

where the second error comes from a t = 1.0 ±0.1, it is indeed not very large. The possible 
theoretical errors induced by the uncertainties of other input parameters are all very small, 
since both z and 5 are stable against the variations of them. 

The theoretical predictions for the direct CP- violating asymmetries A d ^ p (B ± — > p ±r r]^) 
in the PQCD approach are generally larger in size than those obtained by using the QCD 
factorization approach 

A d j P {B ± -> p ± r ] ) = (-2.4 ± 6.4) x 10" 2 , 

A% P {B± - pV) = (4-l^° 9 6 ) x 10- 2 . (79) 

On the experimental side, the new world-average jH| is 

A CP {B ± -> p ± r ] ) exp = (-3 ± 16) x 10- 2 , (80) 

which is still consistent with the predictions in both PQCD and QCD factorization ap- 
proach within the still large experimental error. More data are clearly needed to make a 
reliable judgement. 

We now study the CP-violating asymmetries for B° — > p°r)^ decays. For these neutral 
decay modes, the effects of B° — B° mixing should be considered. For B° meson decays, 
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FIG. 6: The direct CP asymmetries (in percentage) of B + 
(dotted curve) as a function of CKM angle a. 



p + rj (solid curve) and B + 



p + v' 



we know that AT/Am d < 1 and AF/F < 1. The CP- violating asymmetry of B°(B°) 
pOjyW decay is time dependent and can be defined as 



r (g(AQ -> / CP ) - r (gg(At) -» / C p) 
r (lJ(Af) -> / CP ) + r (B°(At) - / C p) 

Agp cos(AmAt) + A™p sin(AmAt), (81) 



where Am is the mass difference between the two B° mass eigenstates, At = tcp — t tag 

is the time difference between the tagged B° (B ) and the accompanying B (B°) with 
opposite b flavor decaying to the final CP-eigenstate fcp at the time tcp- The direct and 
mixing induced CP-violating asymmetries A^ P and A™p can be written as 



-V 7< 2 1 Amix 2Im(\ CP , 



Acp " 1 + |Acp| 2 ' CP " 1 + IW ( } 

where the CP- violating parameter \ C p is 

^(pV)|/r e// |B°) 1 + 

ACP Wa<A (,) l#e//|£ > 1 + **(*+«>■ ^ 

Here the ratio z and the strong phase 5 have been defined previously. In PQCD approach, 
since both z and 5 are calculable, it is easy to find the numerical values of A^p and A™p 
for the considered decay processes. 

For B° — > p°r] and p°r)' decays, the numerical values of the ratio z and the strong phase 
5 are 

z (p°ri) = 4.0, 5(p°r]) = -57°, (84) 
z(p°rf) = 6.8, 8(p°r)') = -65°. (85) 
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FIG. 7: The direct CP asymmetry Aq P (in percentage) of B° — * p°r] (solid curve) and B° — > (Prf 
(dotted curve) as a function of CKM angle a. 

Unlike the case of B ± — > p ± r]^ decays, we here have z > 1, which means that the "P" 
term is much larger in size than the "T" term for B° — > p°r)^ decays, since the "T" term 
here is a color-suppressed tree. 

In Figs. [7| and |H1 we show the a— dependence of the direct and the mixing-induced CP- 
violating asymmetry Aq P and A™ l p for B° — > p°r] (solid curve) and B° — > p°rj' (dotted 
curve) decays, respectively. For a ~ 100°, one can find numerically that 

Afp{B° -> p°rj) « -41%, A^(B° -> p°rj) « +25%, (86) 
A^(5° -> p°r/) « -27%, -> pV) « +11%. (87) 

They are also large in size. The theoretical errors induced by the uncertainties of input 
parameters are only about 10% because of the cancelation in ratios. If we vary a t in 
the range of 0.9 < at < 1.1, however, the theoretical predictions for the CP- violating 
asymmetries of the penguin-dominated B° — > p°r]^ decays may change significantly 

A d S r P {B° -> p°ri) = [-85%, +24%], A™*(B -> p%) = [-19%, +35%], (88) 
^(£?° -> pV) = [-75%, +13%], j4£]?(B -> pV) = [-9%, +22%]. (89) 

This feature may be interpreted as an indication of the importance of the NLO contribu- 
tions to those penguin dominated decay modes. 

If we integrate the time variable t, we will get the total CP asymmetry for B° — > p°rj^ 
decays, 

Acp = ^— 2 Af P + — (90) 
1 + x z 1 + ar 

where x = Am/T = 0.771 for the B° — B° mixing In FigEJ we show the en- 

dependence of the total CP asymmetry Acp for B° — > p°r] (solid curve) and B° — > p°r]' 
(dotted curve) decay, respectively. For a ~ 100°, the total CP asymmetry is around 
-10% for B° -> (Prjff) decays. 
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FIG. 8: The mixing induced CP asymmetry A™p (in percentage) of B° — ► p°rj (solid curve) and 
B° —* p°7]' (dotted curve) as a function of CKM angle a . 
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FIG. 9: The total CP asymmetry A l °£ (in percentage) of B° 
(dotted curve) as a function of CKM angle a . 
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For B + — > p + r] and p + r]' decays, the large CP-violating asymmetry around —15% could 
be measured in the running B factory experiments since their branching ratios are rather 
large, ~ 10~ 5 . For Br(B° — > p°rj) and p°rj' decays, however, it is very difficult to measure 
their CP-violating asymmetries at current running B factories since their branching ratios 
are very small, ~ 1CF 8 only. It could be measured in the forthcoming LHCb experiment. 



D. Effects of possible gluonic component of rj' 

Up to now, we have not considered the possible contributions to the branching ratios 
and CP-violatin g as ymmetries of B — > prj' decays induced by the possible gluonic com- 
ponent of rj' |2H I29I l3p| . When Z, q i 7^ 0, a decay amplitude M.' will be produced by the 



24 



gluonic component of rf '. Such decay amplitude may interfere constructively or destruc- 
tively with the ones from the qq (q = u, d, s) components of rf, the branching ratios of 
the decays in question may be increased or decreased accordingly. 

Unfortunately, we currently do not know how to calculate this kind of contributions 
reliably. But we can treat it as an theoretical uncertainty. For \M' /M(qq)\ ~ 0.1 — 0.2, 
for example, the resulted uncertainty for the branching ratios as given in Eqs.(|68|) and 
(fTOj) will be around twenty to thirty percent. 

From Eq. (|68|). one can see that the theoretical prediction of Br(B + — > p + rf) in the 
PQCD approach agrees well with the measured value within one standard deviation, which 
means that there is no large room left for the contribution due to the gluonic component 
of rf. We therefore believe that the gluonic admixture of rf should be small, and most 
possibly not as important as expected before. 

As for the CP-violating asymmetries of B — > prf decays, the possible contributions of 
the gluonic components of the rf meson are largely cancelled in the ratio. 

V. SUMMARY 

In this paper, we calculate the branching ratios and CP-violating asymmetries of B° — > 
p°r], B° —>■ p°rf, B + — > p + r], and B + — > p + rf decays in the PQCD factorization approach. 

Besides the usual factorizable diagrams, the non-factorizable and annihilation dia- 
grams as shown in Figs. ([!} and (J2J) are also calculated analytically. Although the non- 
factorizable and annihilation contributions are sub-leading for the branching ratios of the 
considered decays, but they are not negligible. Furthermore these diagrams provide the 
necessary strong phase required by a non-zero CP- violating asymmetry for the considered 
decays. 

From our calculations and phenomenological analysis, we found the following results: 

• From analytical calculations, the form factors for B — > rj, B — > rf and B — > p tran- 
sitions can be extracted. The PQCD results for these form factors are Aq~* p (0) = 
0.37, Fq^ v (0) = 0.15 and Fq~* v (0) = 0.14, which are in good agreement with those 
obtained from the QCD sum rule calculations. 

• For the branching ratios of the four considered decay modes, the theoretical predic- 
tions in PQCD approach are 

Br(B + -> p + r] {,) ) ss 9 x 10" 6 , (91) 
Br(B° -> p°rf ,} ) w 5 x 10~ 8 . (92) 

Although the theoretical uncertainties are still large (can reach 50%), the leading 
PQCD predictions for the branching ratios agree well with the measured values 
or currently available experimental upper limits, and are also consistent with the 
results obtained by employing the QCD factorization approach. 

• For the CP-violating asymmetries, the theoretical predictions in PQCD approach 
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are 



13%, (93) 

18%, (94) 

41%, A%£(B -> p°v) ~ +25%, (95) 

27%, A™(B° -> pV) « (96) 

for ct ~ 100°. For .E^ — > decay, the CP-violating asymmetry around —15% 
could be measured in the running B factory experiments. For the neutral decays, 
their CP-violating asymmetries may be measured in the forthcoming LHCb exper- 
iments. 

• The major theoretical errors of the computed observables are induced by the un- 
certainties of the hard energy scale t/s, the parameters Ub and mj, as well as the 
CKM angle a. 

• From the good consistency of the PQCD prediction of Br(B + — > p + i]') with the 
measured value, we believe that the gluonic admixture of rj' should be small, and 
most possibly not as important as expected before. 



Af P (B° -> A) « - 
^p(5°- AO « -- 
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APPENDIX A: RELATED FUNCTIONS 

We show here the function h^s, coming from the Fourier transformations of H^°\ 

h e (x 1 ,x 3 , b ± , b 3 ) = K {^xix 3 m B bi) [0(bi - b 3 )K (^/x^m B bi) I {^/x^m B b 3 ) 

+9(b 3 - bi)K {^/x^m B b 3 ) I {^/x^m B bi)} S t (x 3 ), (Al) 

h a (x 2 , x 3 , b 2 , b 3 ) = K Q (iy/x 2 x 3 m B b 3 ) [9(b 3 - b 2 )K (i^/x 3 'm B b 3 ) I {i^/x^m B b 2 ) 

+9(b 2 - b 3 )K (iy/x^msh) I {i^/xlm B b 3 )} S t (x 3 ), (A2) 

h f (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) = j#(& 2 - bi)l (M B y/x 1 x 3 b 1 )K (M B y/x 1 x 3 b 2 ) 

> ( K (M B F(!)6i), for > 0\ 

+ (^-M)^ fH a) (MB ^ 6iUor Fh<0 j, (A3) 
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hf(xi,X2,X3,bx,b 2 ) = - b2)Ko(i^X2X 3 b 1 M B )l (iy/x2X 3 b 2 M B ) 

} f K (MbF( 2 )6i), for F ( 2 2) >0 N 
+ ^"^\\fEt\M B ^\b 1 ), for Ff 2) < )' < A4 > 

hj(xi, x 2 , x 3 , 61, 6 2 ) = \ 0(h - b2)Ko(i^X2X 3 b 1 M B )lo(i^X2X 3 b 2 M B ) + (61 <-+ b 2 ) 



7T«, 



Hp (V^i + x 2 + x 3 - - x 2 x 3 61 Mb) 



(A5) 



where Jo is the Bessel function and Kq, Iq are modified Bessel functions Kq{— ix) 
— (n/2)Y (x) + 1(71/2) J (x), and Fq^s are defined by 



F {2) 



(xi - x 2 )x 3 , 
(xi - x 2 )x 3 . 

The threshold resummation form factor St(xi) is adopted from Ref.|2 

2 1+2 T(3/2 + c 



S t x) 



V5FT(1 + c) 



[x(l-x)] f 



where the parameter c = 0.3. This function is normalized to unity. 
The Sudakov factors used in the text are defined as 



S a b(t) = s (xxm B / V2, bi^j + s (x 3 m B / y/2, b 3 J +s((l- x 3 )m B / V2, b 



1 



In 



ln(t/A) 



In 



ln(t/A) 



S cd (t) 



ln(6iA) ' "~-ln(6 3 A) 
s [x 1 m B j y/2, bA + s (x 2 m B / \/2, 6 2 J + sUl — x 2 )m B / y/2, b- t 

+s (^x 3 m B j\pl^b\ 

) + s ({1 - x 3 )m B /V2,b^j 

1 



Pi 



2 In ln(t/A) + ln ln(t/A) 
-ln(6iA) -ln(6 2 A) 



Sef(t) 



s X\m B 



2, bij + s (x 2 m B / V2, b 2 J + s f (1 — i 2 )m B / a/2, 6s 
+s (x 3 m B / V2, b 2 J +s ({1- x 3 )m B / a/2, 6 5 



1 

"A 



In 



ln(t/A) 



2 In 



ln(t/A) 



ln(6iA) -ln(6 2 A) 
= s (x 2 m B / y/2, 61 J + s {x 3 m B / V2, b 2 J + s ^(1 — x 2 )m B / \/2, 61 



ln ln(t/A) + ^ ln(t/A) 



(A6) 
(A7) 



(A8) 



(A9) 



(A10) 



(All) 



(A12) 



ln(6iA) -ln(6 2 A) 

where the function s(g, 6) are defined in the Appendix B of Ref. Q and the hard energy 
scale t'jS have been given in Eq. (J3l)j) . 
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